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The volume dependence of the ground-state energy of n identical bosons with short-range inter- 
f~>) ■ actions in a periodic spatial volume with sides of length L is calculated at order L~ 7 in the large 

| volume expansion. This result will enable a refined determination of the 7r + 7r + 7r + interaction from 

lattice QCD calculations. 
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It is now well established that two-body interactions between hadrons can be studied with lattice QCD as the 
volume dependence of the energy spectrum of two hadrons is related to their scattering amplitude below inelastic 
thresholds [H, 0. Recently, this method has been used to determine the tt + tt + scattering length a 7r + 7r +, with 
~ 1% precision with a rif = 2 + 1 fully-dynamical mixed-action lattice QCD calculation. In order to extract the 
many-body interactions from lattice QCD calculations, the energy of multi-hadron states in a finite volume must be 
calculated with lattice QCD and combined with the known dependence of this energy on the many-body interactions. 
The ground state energy of a system of n identical bosons with short-range interactions in a cubic volume with sides 
of length L was recently computed at (D(L~ 6 ) in the large-volume expansion The underlying motivation for that 
— — 1 , work, which builds upon the classic works of Refs. H |g, 0, , was to provide a way to determine the three-body 
interactions between 7r + 's from lattice QCD calculations, which first enters at that order In Ref. @, this result 
was used in conjunction with lattice QCD calculations of multi-pion systems to determine the interaction between 
■ three 7r + 's for the first time. In order to refine the determination of the tt + tt + tt + interaction, here we compute the 
\Q | contribution to the energy-shift of n identical bosons at 0(L~ 7 ) in the large volume expansion. The energy-shift of 
three identical bosons in a finite- volume has been been computed recently in Ref. [lfj, and our n — 3 calculation 
' agrees. 

The ground-state energy of n identical bosons is calculated using standard Schrodinger perturbation theory, with 
a Hamiltonian, appropriate to the order we are working in the large volume expansion, of the form 
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where the operator annihilates a tt + with momentum k with unit amplitude. The divergences that arise at loop-level 
are regulated with dimensional regularization, and therefore the coefficients of the two-body interaction can be readily 
identified with the parameters describing the scattering amplitude: the scattering length, a, and the effective range, 
r (pcotJ = — i + \rp 2 + ... ). The terms proportional to M~ 3 in eq. (|TJ) describe the leading effects of relativity. 
Only the momentum independent three-body interaction, %(/^), is required at 0{L~ 7 ). Our method of computation is 
equivalent to the pionless EFT describing low-energy nucleon-nucleon interactions, EFT(7f) [ll|, [Tl, EH (when modified 
to describe systems with natural scattering lengths) and the method of pseudo-potentials used in our previous work. 
The divergences that occur in loop diagrams are renormalized order-by-order in the expansion, preserving the power 
counting, and hence the explicit dependence of the bare three-body coefficient on the renormalization scale, /i. 

The calculation of the energy-shift of n identical bosons at 0(L~ 7 ) due to the interactions defined in cq. ([T]) is 
straightforward but tedious. We will not delve into the details, referring the reader to our previous work [H and that 
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of Ref. [l0|, and simply state the result. The energy-shift of the ground state is 
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where the geometric constants that enter are 1 
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and (fc)=n!/(n — fc)!/fc!. The last term in the last bracket of eq. (|2|) is the leading relativistic contribution to the 
energy-shift. Deviations from the energy-shift of n-bosons computed with non-relativistic quantum mechanics arise 
only for three or more particles as the two-particle energy-shift has the same form when computed in non-relativistic 
quantum mechanics and in quantum field theory [l],^]. I n ec L- ©j ^ms is the value of the scheme-dependent quantity 
S in the Minimal Subtraction (MS) scheme that we have employed to renormalize the theory (a change in scheme 
results in a change in S and a compensating change in 773 (/it V). 2 The % are renormalization scheme independent. Our 
result at n = 2 agrees with large volume expansion of Ref. [y, [1J , and at n = 3 agrees with the previous computation 
by Shina Tan [Tof. 

The renormalization-scale independent, but volume dependent, quantity 
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was determined in recent lattice QCD calculations 0. It was found to be non- vanishing in systems of three, four 
and five 7r + 's at pion masses of ~ 290 and 350 MeV in a ~ (2.5 fm) 3 volume, when extracted at 0(L~ 6 ) in the 



1 The constants X, J , K, were defined previously in Ref. f3|, while the constant C is defined to be the integer-triplet sum 

n^O 1 1 

and is equal to C = cm in the notation of Ref. llOH . The constants T^i arise from combinations of up to three-loop diagrams, and involve 
three-, six- and nine-dimensional sums over integers, and can be written in terms of constants defined in Ref [10| plus one additional 
sum, Si, 

To + Tin = ^aiAAi-laiAi + -(2n-9)a 2 Ai + ^(n-4)aiBi--(7n-29)£ + 2(n-3)Si , 



where 



\ ^ = qo 4221 Fi 
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In the notation of Ref. [J, 5ms = 2<2 + 7S. The numerical value in eq. J3J corrects a minor error in Q in Ref. Q|. 
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large volume expansion. Its size was found to be consistent with expectations based upon naive dimensional analysis, 
773 ~ l/(m,r/^). Our result will allow for further refinement of such extractions. 
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